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CALCUU^ION OF TRAHSONIC PLOWS USIHG AH EXTEHDBD INTEGRAL EQUATION 
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Abstract 

An extended integral equation method for transonic flows is 
developed. In the extended integral equation method velocities in 
the flow field are calculated in addition to values on the aerofoil 
surface, in contrast with the less accurate ’standard* integral 
equation method in which only surface velocities are calculated. 

The results obtained for aerofoils in subcritical flow and in 
supercritical flow when shock waves are present compare satisfactorily 
with the results of recent finite difference methods. 


* Senior Research Fellow 
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lairoduotion 

At the present time the most common method of calculatiiMl the 
inviscid pressure distribution around aerofoils at transonic speeds 
is by using finite difference techniq.ues. For subcritical flows 
one of the most accurate methods is that of Sells (l). For super- 
critical flows with shock waves the initial work of Murman and Cole 
(2) has been followed by many developments, notably those by Kuriaiai 
and Krupp (3), Garabedian and Korn (4) and by Jameson (5). It is 
pointed out by Murman ( 6 ) , however, that the difference scheme 
used in these methods is incorrect at the shock location, leading 
to incorrect shock jump relations; corrected results are presented 
in Ref, ( 6 ), 

One of the earliest attempts to solve the transonic flow problem 
was the integral equation method first developed for non-lifting 
flows by Oswatitsch (7) with subsequent extensions by Spreiter and 
Alksne ( 8 ). Later developments to include lifting flows are by 
Norstrud (9) and Nixon and Hancock (lO) although it should be noted 
that the formulation used by Norstrud for lifting flows is incorrect 
(ii). 

In the integral equation method for tv;o-dimensional flows the 
non-linear partial differential equation for the perturbation vel- 
ocity potential is written in integral form using Green* s Theorem, 

The resulting integral equation for the velocity potential can 
be differentiated to give the perturbation velocity in terms of 
a line integral over the aerofoil chord, involving only linear terms 
and a surface integral over the flow field involving only non-linear 
terms. The line integral can be easily evaluated using standard 
methods but the evaluation of the field integral requires a know- 
ledge of the velocity distribution over the Entire flow field. In 
the standard integral equation method the field integral is reduced 
to a line integral by representing the transverse variation of the 
velocity by a fairly arbitrary approximation function involving 
only the velocity on the aerofoil surface. The resulting integral 
equation for the surface velocity can he solved without much diff- 
iculty, This approximate evaluation of the field integral, although 
accurate for shock’ flows, is not sufficiently accurate to 
give satisfactoiy resmlts when shock waves are present. 


Many problems in subsonic * linear* aerodynamics, both in stea^iy 
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flow oscillatory flows, oaa be formulated in terms of integral 
equations and there exists a considerable expertise in tbe solution 
of these equations. Similarly many problems in transonic aerodyn^ics, 
again both for stea^ and oscillatory flows, can be formulated as , 
integral equations! examples of steady flow problems are given ly 
Norstrud (9) and an example of oscillatiiy flows is given by Nixon 
(12), fhe trsmsonic integral equation contains a line integral 
similar to one appearing in subsonic aerofoil theory and the available 
expertise in the evaluation of this integral can be incorporated into 
susy solution of the transonic integral equation, fhe main advantage 
of the integral equation methods over finite difference methods is 
that, in principle, the numerical solution of the integral equation 
is easier than the solution of the differential equation and is 
typified by rapid convergence. The main disadvantage of existing 
solutions of the transonic integral equation is the unsatisfactory 
results due to the rather primitive and inaccurate approximation 
of the field integral component of hhe transonic integral equation. 

If the accuracy of the evaluation of the field integral can be 
improved then it is suggested that integral equation methods will 
become a useful technique in a wide range of problems in transonic 
aerodynamics. 

In the extended integral equation method presented in this paper 
an alternative means of evaluating the field integral is developed. 

The flow field is divided into a number of streamwise strips and the 
transverse variation of the perturbation velocities across each of 
these strips is approximated by an interpolation function in terms 
of values on the strip edges. The field integral is then reduced 
to a line integral which is in turn evaluated by quadrature. The 
fundamental integral equation is thus approximated by a set of non- 
linear algebraic equations. 

The pressure distribution around both lifting and non-lifting thick 
aerofoils in subcritioal flow is calculated and compares favourably 
with the results of o^her exact methods. The pressure distribution 
around aerofoils in supercritical flow with shock waves is also 
calculated and there is good agreement with the results of recent 
finite difference methods. The computing time for a lifting sup- 
ercritical example is about 8O-IOO secs, on an ICL I9O4S computer. 
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B^io Bquatlons 

A two dimensiotial cartesian co-ordinate system is chosen with 
the origin at the leading edge of the aerofoilf the x-azis is in 
the freestreaffi direction ^d the s-^is is normal to the freestres^l 
X and m are non dimensionalized with respect to the aerofoil ehord« 

In order to sirapli:^ the govenaing eq^aations for transonic flow 
it is frequently assumed that the flow is hoth isentropic and irrot- 
ational even when shock waves are present, Thj perturbed flow 
may thus he represented hy a perturbation velocity potential, 

X, z), defined by 

« H- (» I 

where u and w are the non-dimensi 4 nal perturbation velocities in 
the X and z directions respectively, relative to the freestreara 
velocity. The full potential equation for steady flow is of third 
order in ^ but for transonic flows a number of the higher order 
terms may be neglected for post practical applications. 

The variables 

I y * X j z-/»z (2) i 

are introduced where, if is the freestream Mach number and Y 
is the ratio of specific heats 

A = 

and 

k » (3) I 

A second order potential equation for transonic flow can then be 
written as 

where r 

I 
} 

I 

If the second term on the right-hand side of Eq, (5) is neglected 
then Eq, (4) reduces to the transonic small disturbance equation. 


= (♦) 
s (s) 



The boundary conditions arei 
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(a) that tha parturhation valocity potential and its derivativas 
vanish at infinita distance upstreaa of the aerofoil, 

(h) that the flow direction of the aerofoil surface is tangential 
to the aerofoil surface, 

(c) that for a suhsonic trailing edge the Kutta condition of finite 
pressure at the trailing edge must he satisfied* 


If z ss Zy(x) and z = z^(x) are the equations of the upper and 
lower surfaces of the aerofoil respectively then the tangency 
boundary condition can he writt en in the vari ahles of E g* (2) as 


w(S.Ja= Z^V[n-f'ikx,zjf 
wa,^)= (6j 

where 

2,W= 


It can he shown tb'^t Eg, (4) is elliptic (suhsonic) when 
hyperholic (supersonic) when O^z-C^lpS^onic conditions exist when 
o=/-(vr When Eg, (4) is hyperholic in character then a solution 
of Eg. ( 4 ) may have discontinuities in the velocity components u,w. 
These discontinuities are the shock waves of Eg, (4). The shock 
jump relations can he found fairly easily from Eg. (4) since it is 
writer en in conservation form* thus 


f ^ ■ 9 V + [ = 0 (7«) 

where ^ J^denotes the jump across the shock wave and cos (n,^), 
cos (n,^ ) are the direction cosines of the shock wave. It should 
also he noted that 



(7b) 


Treatment of Boimdary Conditions 


There are two common methods of treating the houadary conditions 
in aerofoil theory, namely thin aerofoil theory and thick aerofoil 
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■feheoiy. Thin aerofoil theory introduces an approximation to the 
exact prohlem using the assumption that the disturbano 'S from 
the freestream are of the same order of magnitude as some power 
of the thickness/ chord ratio of the aerofoil section, which is 
usually small compared with unity. The boundary conditiens and 
aad the potential equation are then expanded as a series in powers 
of the thickness/ohord ratio giving an infinite set of differential 
equations with their associated boundary conditions . Generally 
only the first few terms of these series are retained. The tangenoy 

boundary condition for each equation Is now satisfied on the plane 
z = 1 0 rather than on the aerofoil surface. The accuracy of thin 
aerofoil theory can be improved by progressively including more 
terms in the series. This method of solution gives rise to sing- 
ularities due to the assumption of small disturbances breaking down 
in the neighbourhood of the leading edge and these singularities 
ira;;t be removed in the course of the solution. 

In t ick aerofoil theory no assumption regarding the magnitude 
of the flow disturbances is made and consequently thick aerofoil 
theory is generally more difficult than thin aerofoil theory and 
the numerical procedures are more complex although, of course, the 
•exact* result is computed. 

The considerable experience available in using both the’steindard* 
(8,9,10) extended integral equation method is based on the 

formulation of the problem using the thin aerofoil boundary conditions, 
particularly that the boundary conditions are satisfied on the plane 
z =±0 rather than on the actual aerofoil surface. In transonic 
flows when the non-linear potential equation, Eq,(4)» is applicable, 
the series solution normaly used when the non-linear terms in Eq,(4) 
can be neglected is impractical and some alternative must be sought. 

The principle requirement is that the tangency boundary condition 
must be satisfied on the plane z 

If it is assumed that the external flow (that is, external to 
the aerofoil) can be continued analytically inside the aerofoil 
to the place z - +0, z « -0, for the upper and lower half planes 
respectively, with the corresponding boundary conditions on the 
plane z =^0, then existing expertise can be used in formulating 
the integral equation for a modified problem, the domain of which 
now includes both the real external flow and a fictitious internal 
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flow. Since It is assumed that the interior flow is an anal^io- 
continuation of the external flow then the value of can 

he obtained from the physical hounda^T- conditions B<i, (6) using 
a Taylors series expa nsion ! thus for example is given hy 

a (8J ! 

If only the first two terras of the series in Iq,. (8) are taken then 
on using Sqs, (4f6) the boundary conditions for the modified flow 

s2jp^*f(ktxH 

( 9 ) 

The boundary conditions, Sq. (9)» used in the subsequent 
analysis. Since there is no assumption of small disturbances 
except possibly in the number of terms retained in the Taylors 
series expansion, no singularities due to approximations in the 
boundary conditions need arise. Once the solution to the modified 
problem is known the solution to the physical problem is found by 
considering only external flow. 


The Integral Equations 


The partial differential equation for the velocity potential 
Sq, ( 4 ), with the associated shock relations Eq, (7)» can be writtem 
in integral form using Green* s Theorem! thus 

where ( r, f* ) are co-ordinates equivalent to(***^ , 

and the operator ’ A * is defined for as 

AfCt)s (12) 

The surface integral ih Bq.tlO; is defined J'or 3T>0 as 

where 3^ , denote the location of the shock waves on the upper 
and lower half planes respectively. 






Iq. (lO) CM be differentiated with respect to M to give the 
velocity component thus 




The term 


gfr li) 


+aj|£ - 

^ (14) ' 

I 

arises from differentiation of the limits around 


the singular point in the field integral. If the second 

integral on the right-hMd side of Iq, (I4) is integrated by parts 
then , 




Qi-) 


It may be shown that 


Inhere 


S' 06) 

I/8.S.S,) 


(W)i 


and the surface integral in Eq, (17) is_ defined by 

rW ,te 




On substition of Eqs. tl5f 16) into Eq, 

Oca,«- 8c».« = l^a.z) -f X,ca,s. V) 


(19) 


where 




( 20 ) 


If the limit as S**+0 is taken it can be shown that Eq, (j.9) 


reduces to 


{0CS,4C»- 


which is symmetric with respect to the 2— axis, A similar result 
occurs if the limit as ^■♦—0 is talten. The additional equation 
to give the anti-symmetric component of UC^^to) , A iJ(i) , can be 
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oMained 'by differeatiating Eq, (lO) with respect to 1 and taking 
the limit as Thus , after some manipulation 


+XjLi.V , 22 ) 


where 




(23) 


and 


where 


S (24) 


?n 

Eq. ( 22 ) can bs inverted to give 


(2S) 


Combination of Eq.(25) and Eq,(2l) then gives the required 


ie redefined as 


If the integral Il(-^ ±0) 

!.«.*«= J*PWS3SW)*% 

and if the integral is redefined as 


(27) 


then Eq. (l9) can be uedd for all -f . 


In order to simplify the present analysis it is assumed that 
any shock waves in the flow are normal to the freestream. Using 
this assumption the shock relations, Eq, (7a) become 
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If shock waves are present in the flow then, with g(*i*) given 
by Eq. (5), it can be shown that in order to ensure finite accel- 
eration through the line (0 ^ l), Eq, (l9) must be solved 

subject to the conditions 


1 


*x.»,«) * 2 ^ 0 . 1 , Mil 




M 

X 




•o 




(^9a) 

un) 


where denotes the iO^tJ line. It c«i he shown that Eq,. (29c) 

is alweyg aatisfledf The remaining equations Eq,(29a, 29h) are 
sufficient to give location of the line X«Ci> aiid the shock 

loc&tione Hpi/f ^4 • 


The perturhation velocity component UCXf2) oan then found 
from Eq.(jl9) with Eq,(9) &nd Eq,('d9). Having found U<^g2J the 
velocity component can he found using the irrotationality 

relation 

ao) 


C£'CS,X)= ucr^tidf 


Evaluation of the Int e grals 

In order to solve Eq,(l9) the line integral and the 

field integrala in need to he evaluated. The line 

integral involves only the value of on the aerofoil surface 

through the hourdary conditions Bq,(9) and can he evaluated using 
standard techniques in terms of specific values of UL9^'i) on the 
aerofoil surface. The evaluation of the field integrals in 
X.C9,*,St,) is more cc ,«plex since they involve 3c«f) over the field 
An approjdmate, although sufficiently accurate, evaluation of 
) is possihle if 9^4 f) is known at specific points in 
the flow field and an interpolation function used to express 
over the rest of the flow field. The practical requirement is an 
accurate estimation of the velocity on the aerofoil surface and hence 
an acceptahle solution to Eq.(iy) neid not necessarily be very accurate 
in the far field provided the surface velocities are calculated to 
sufficient accuracy. 

Let the flow field he divided into 2!f strips as shown in Fig*l, 

In principle the outeimjost strips in each half pl^e extei^ to infinity 
hut in practice they need only extend to a stitahle finite location 
in the. far field. It is assumed that the function is 

known on each strip edge i fs ^ ^ ~ ZA/ ) and ^^e variation of 


10 


$if 9 ^ ) is tach itrip is rsprsesntsd ’by’ sons intsrp^la'^ion flmctioni 
in tne present work linear interpolation is used. Thus in the 
strip 

I jcgpK >-f,) 

llN / 

5 


Using Sq,, (3l) the integrations with respect to f in the field 
integrals ciai he perforoed. The suhsequent line integrals in the 
f direction can he evaluated using standard methods in terms of 
si^if ) where the (i « 1»^) are specified values ofw. On 
putting 1? and * ■ the integral equation Iq, (19) reduces 

to a set of non-linear algebraic equations for U , ^ ) if 
^ ) is assumed kno^f thus 


( 32 ) 


which ma^ he solved "by iteration. 


je the 0 f are known the (Sj ) can he 

(30) and the surface velocit^eB can he foun'1 using 


found frorri 
interpolation. 


^ince high accuracy is required only at the aerofoil surface 
the function on e. certain number of outer strip edgec 

may be adequately estimated in terras of values on the inner strip 
edges using some suitable approximation function without eignific- 
antly affecting the accuracy of the calculation of the surface 
velocities* In the present work it is assumed that 



wnere 




and ? is some value of ? in each half plane. 


Calculation Procedu re for Sup ercritical Flow with Shock Havea 

When no shook waves are present Eq, (32) can he e-olved directly. 
When shock waves are present Iq, (32) must he solved subject to the 
conditions lq,(29). It is difficult to sotlsfy Iq, (29) at each 
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stf^e of the iter&tioa by alterl^ the shook locetioa since the 
location of the discontinuity in oCSIyl) then changes at each 
iteration* An alteraati've aethod of satisfying Iq., ( 29 ) is as 
follows. 


E«l.(32) can he tfrit ten in th e form_ _ 


where the S (l|^ is a par^eter which is constant along each strip 
edge. An initial estimate is ma4e of the shock locations , 

Xn f the function ) and 1?he integrals 

5^) and 2^ , 5^ ) can then he evaluated, fhe parametsr 

f (if) and the locat ion of thefO= 0 point on the strip edge 
(2T= ^e evaluated hy enforcing the condition of finite 

streamwise acceleration along each strip edge, that is satisfying 
Eq, ( 29 a, 29h) along each strip edge. New values of U (x,* , ), 

si^i , ^i) are then computed and the prscess repeated until conver- 


gence of the f(^). The shock locations are then changed and 
the procedure outlined above repeated. The shock locations are 
correct when the £ (2,*) are ixnity. Since high accuracy is re- 
quired only on the aerofoil surface the accuracy of the solution 
in the flow field can decrease as /r/ -^ao provided the surface 
solution is not affected to any significant degree. In practice 
this is taker '.o mean that the converged value of £(2^) may he 
allowed to differ from unity on the outermost strip edges. 


Results 


The pressure distributions over several aerofoils have been 
computed for both subcritical and supercritical flows with shock 
waves. For the subcritical flows the function is defined 

by Eq, (§) and the transonic parameter /?, Eq,(3)> is defined after 
Hancock to be 

(3r) 

For supercritical flows with shock waves the function g(K,S) 
is taken to be 

in order to siraplify the computation,, and the is 

taken to be s, 

Of <57; 


In all caaes the exact pressure-velocity relation is used. 
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Three es:amples of suhcrltloal flo^ caleulations are shown in Figs. 2-^ 

and are compared to exact results. The results of the standard integral 
fl5) 

equation method ^ are shown for comparison. In Fig, 2 the flow around 
a NACA 0012 aerofoil at 0^ incidence luid » 0.72 is shown smd it can he 
seen that the present results compstre well with the * exact* results of 
Sells In Fig, 3 the pressure distribution around a HIiR 0,1 - 0,75 

- 1,25 aerofoil at 0® inoidenoe lynd » 0.745 is shown. The agreement 
'.:ith the * exact* result is fairly good except in the vicinity of the 

leading edge. In Pig, 4. the pressure distribution around a HAGA 0012 
aerofoil at 2® incidence and » 0,63 is shown. The agreement with the 

(16) 

* exact* result'* * is good again except in the neighbourhood of the 
leading edge, A similar agreement with the ‘exact* result is shown in 
Pig. 5 * f’or the pressure distril ition around a 14^ th- k NPL 3111 aerofoil 
at 1,2° incidence ana M„ = 0.66'|, 

Three examples of supercritical flow are shown in Pigs, 6-8. In Fig, 6 
the pressure distribution arouid a NACA ool2 aerofoil at zero incidence 
and 0 ,816 is shown and the present results agree fairly well with 

the finite difference results calculated using the method of Garabedian 
and Korn^"^^ us’ % the non-conservative difference scheme. In Fig. 7 the 
pressure distr \ition around a HACA 0012 aerofoil 2 ° incidence and M^=0,75 

(17) 

is compared to the finite difference results calculated Bauei* and Korn' 

using both ‘non-conservative* and ‘conservative* difference schemes. As in 
the non-lifting example the present results agree fairly well with the ‘non- 
conservative* results. In Pig. 8 the pressure distribution around a NACA 
64A4IO aerofoil at 0° incidence euid M*, = 0,72 is shown. As in the previous 
example the present resul^ agree fairly well viith the ’non-conservative* 
finite difference results. This tendency of the present calculations to 
agree with the * non-cons ervstive* rather than the correct ‘conservative* 
results may be attributed to the use of the transonic small perturbation 
ecuation for the integral eq,uation calculations as opposed to the use of 

(17^ AX (1®) 

the full potential equation by Bauer and Korn ' ,ana Jameson » 


For the subcritical examples 30 streamwise elements and 5 strips in each 
half plane are used' , The velocity on the outer strip edges in each Iialf 

plane is estimated using Eq,(33), Convergence is obtained after about 10 
iterations. 

In the supercritical examples 30 streamwise elements and 8 strips in each 
half plane are used. The ve’^ocity on the 3 outer strip edges is 

estimated using Eq.(^j). ?<-i each shock location 


convergence of the €{^) is obtained after about 10 iterations, 

A tn>ioal computing time for the lifting supercritical ex^ple, 
is 80 - 100 secs, on an ICL 1904S, Further improvements in the 
numerical scheme are feasible and it is probable that the coraputii^ 
time can be considerably decreased, 

Gonolusions 

An ertended integral equation method has been developed for the 
transonic flow around lifting euid non-lifting aerofoils, The 
pumber of iterative steps to convergence is small and this is 
reflected in the comparatively low computing time. Improvements 
in the numerical solution of the integral equation arc possible 
which should decrease the overall computing time. The results 
obtained for both lifting and non-lifting aerofoils in subcritical 
flow agree satisfactorily with ’exact* results. The extended 
integral equation method gives considerably improved results over 
the earlier standard integral equation method. The results for 
supercritical flows with shock waves are in fairly good agreement 
with existing finite difference results, although there appears 
to be better agreement with the incorrect ’non-conservative* 
finite difference results than with the correct ’conservative’ 
results. This is attributed to the present use of the trsmsonic 
small disturbance potential equation as opposed +o the full 
potential equation used in the finite difference calculations. 
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Arrangement of strips. 
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Pressure distribution around a NACA 0012 aerofoil at 
0® Incidence; * 0.72 . 

Fig. 2 



Pressure distribution around an NLR 0.1 - 0.75 - 1.25 
aerofoil at 0® Incidence; Moo - 0.745 

Fig. 3 



Pressure distribution around a NACA 0012 aerofoil at 
2® incidence; M^ = 0.63 

Fig. 4 
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Prtssure distribution around a NACA 0012 aerofoil at 
2° Incidence; * 0.75 

Fig. 7 
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Pressure distribution around a NACA 64A410 aerofoil 
at 0° Incidence; **0.72 

Fig. 8 


